How Ice enables Superconductivity in Na x Co02-yH 2 by melting charge order: 

Possibility of novel Electric Field Effects 
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Charge ordering in doped C0O2 planes near the commensurate fillings x = \ and | are con- 
sidered for Na x CoC>2-yH20 and suggested to be competitors to superconductivity, leading to the 
experimentally seen narrow superconducting dome bounded by commensurate doping: 1 < x < |. 
Intercalated hydrogen bonded H2O network, by its enhanced dielectric constant, screen and frustrate 
local charge order condensation energy and replace a generic 'charge glass order' by superconductiv- 
ity in the dome. An access to superconductivity and charge order, available through the new water 
channel, is used to predict novel effects such as 'Electrical Modulation of Superconductivity' and 
'Electroresistance Effect'. 



Discovery of superconductivity in Na x CoC>2-yH20 by 
Takada and collaborators [1] have opened the possibility 
of realizing unconventional superconductivity and novel 
quantum states in 2D arising from strong electron corre- 
lations in doped C0O2 layers. Water of a right proportion 
(y w I ) seems absolutely necessary [2-5] for stabilizing 
superconductivity, suggesting f^O's critical role. While 
water does wonders in nature, its key role here is some 
what puzzling. Elucidating its role in this unusual su- 
perconductor is an important task from material science 
and physics point of view. This is what the present paper 
attempts using phcnomcnological and theoretical consid- 
erations. 

Enthused by the remarkable discovery of superconduc- 
tivity in Na x Co02-yH 2 , the present author [6] and 
others [7] have suggested a single band t-J model as an 
appropriate model to understand superconductivity and 
low energy electronic phenomena. A phase diagram has 
been suggested using ideas of resonating valence bond 
(RVB) theory developed for cuprates [8]. A recent ex- 
periment [3] which shows superconductivity in a rather 
narrow range of doping 7 < x < -|, than predicted by 
RVB theories, suggest that there are perhaps left out in- 
teractions and consequent competing phases which make 
a simple t-J model valid only for limited range of x. The 
situation is not unusual - even in cuprates a simple t-J 
modeling is strictly valid only in the neighborhood of op- 
timal doping. Charge order phenomenon [9,10] is known 
in cuprate superconductors; and it has been suggested to 
compete [11] with superconductivity. 

Doped C0O2 , compared to CuC>2 layers of high 
T c cuprates, has a narrower conduction band and a less 
polarizable valence band of oxygen. Consequently, short 
range coulomb repulsions among carriers are screened 
less. This is likely to stabilize a variety of frustrated 
charge ordering in the triangular lattice, as we discuss in 
this paper. NMR result of Ray et al. [12] indeed provides 
a first evidence for charge freezing in Na x Co02 family 
(x = |), Nao.sCo02 below about T « 300iv, a large 
temperature scale. 



We estimate and include unscreened short range 
coulomb interactions in the t-J model for the study of 
unhydrated Na x Co02 ■ We show that in addition to su- 
perconductivity, charge ordering in the narrow conduc- 
tion band of C0O2 layer is a major instability, for a range 
of higher doping than suspected (Wen et al. in ref 7) . As 
the unscreened near neighbor coulomb interactions are 
large and comparable to the band width, the character- 
istic charge order temperatures are T c h ~ 400JT. 

Fortunately, H2O , in hydrated Na x CoC>2-yH20 , 
makes t-J modeling valid for a range of doping. For rea- 
sons which we elaborate in the present paper, hydrogen 
bonded H2O dipoles of the ice layers screen and frustrate 
charge order condensation energy. That is, they effec- 
tively screen out short range repulsions, and enable the 
physics of a simple t-J model to be realized in a narrow 
range of doping, j < x < 3, as superconductivity. 

We discuss few important charge ordered states at 
commensurate fillings, x = \ and |, which we believe are 
competitors to the experimentally observed superconduc- 
tivity, in the range | < x < |. These reference charge 
ordered states are strongly frustrated by the random po- 
tential from the neighboring Na layers, resulting in a 
glassy phase in the region j < x < ^ and beyond. The 
charge glass phase is likely to be an anomalous metal, 
very much like the spin gap phase in cuprates, where 
there are local charge order activities at low frequency 
scales. 

We estimate the enhancement of the background static 
dielectric constant at short distance due to hydrogen 
bonding in the H 2 layers. We find that this screening 
is sufficient to reduce the large charge order transition 
temperature down to ~ IK and allow superconductivity 
to emerge. Strong commensurability effects and the asso- 
ciated short range charge order reduce superconducting 
T c considerably as we approach the commensurate ends 
x = j and |. 

As H2O stabilize superconductivity and discourages 
charge order, we have a new access to the electronic 
phases of C0O2 layer through water. This leads to the 
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possibility of some novel effects: i) 'Electrical Modula- 
tion of Superconductivity' by external electric field or 
microwave radiation and ii) 'Electroresistance Effect' in 
the normal state. We estimate that voltages <~ 500U, 
applied capacitively to thin films of Na x Co02-yH 2 of 
thickness ~ 1 micron will orient the water dipoles and 
reduce the short distance dielectric screening, resulting 
in stabilization of charge glass order phase and destabi- 
lization of superconductivity. This interesting switching 
effect may have device potential. 

Recently we modeled the low energy physics of doped 
C0O2 using a t-J model and discussed an RVB scenario 
for superconductivity including a PT violating d\ + id 2 
wave superconductivity and a pi + ip2 wave supercon- 
ductivity at a higher doping. To study charge order we 
must include some leading short distance carrier-carrier 
and carrier- iVa-ion screened coulomb interaction: 

Huv = ~tJ2cLC ja + H.c. + JJ2& ■ Si - \rnnj) 
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Here C's and S's are the electron and spin operators. 
As we have an electron doped system we have the 'zero 
occupancy' constraint ^2 a n ia 7^ at every site i. 

Recall that doubly occupied Co 3+ sites carry a charge 
— e with reference to the neutral C0O2 layer and is 
the screened coulomb repulsion between them. We have 
ignored the small two body off-diagonal coulomb interac- 
tion terms. For practical purposes only the nearest and 

next nearest neighbor terms V\ « — % — e A »b and V 2 



are important. 



Here e„ 



£0 + £h 2 o 



represents the short distance dielectric screening arising 
from the filled oxygen bands of C0O2 layers and H2O 
layers in Na x Co02-yH 2 . And \ a b w Co — Co distance 
is the Thomas Fermi screening length for our tight bind- 
ing metallic layer. The random site energy et of charge 
degree of freedom represents the screened coulomb at- 
traction from neighboring Na + ions. 

Electronic structure calculations [13] give a value of 
t w —0.1 eV for the conduction band of the C0O2 layer. 
We estimate V\ and V 2 for Na x CoC>2 , the non-hydrated 
case. The dielectric constants of oxides of Fe and Ni that 
flank Co in the periodic table are ~ 4 to 12. We assume 
a background (short distance) static dielectric constant 
of e a b ~ 8 for our C0O2 layer. Recall that in cuprates 
the background e is large ~ 20, in view of a more polar- 
izable octahedral oxygen network; in C0O2 the oxygen 
filled band is less polarizable and relatively deep below 
the fermi level. Using this dielectric constant and values 
of Co-Co distances in Na , 5 CoO2 we get V\ « 0.8 eV and 
V 2 ~ 0.4 eV. The mean square fluctuation of the carrier 
site energy due to disordered iVa-ions is {5ef) « 0.2 eV 

In the absence of hopping, the dopant carriers Co 3+ 



will order classically and undergo order-disorder tran- 
sition at a fairly high temperature fcsT c h (classical) w 
2V - 10 3 K; here V = \{Vx + V 2 ) is a mean short dis- 
tance repulsion. However the electron dynamics reduce 
T c h (classical) considerably To estimate this reduction 
we perform a mean field analysis of the t-J-V model for a 
CDW order, pretending that an unfrustrated charge or- 
der arises from nesting instability for Nao.$Co02 ■ This 
gives us a BCS like expression for T c : 
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Here p is a fermi sea averaged particle-hole density of 
states corresponding to the ordering wave vector. Sub- 
stituting ep w 0.5 eV, po w — and V ~ 0.4, we get 
T c h w 480K. Frustration on the triangular lattice at half 
filling and disorder effect from Na ions will further re- 
duce this. Thus we get a charge order temperature in 
the right range, T ch (NMR) w 300 K, seen in NMR, 

x= 1/4 




FIG. 1. 'Classical' charge order for x = |. Spin-0, charge 
— e carriers (Co 3+ ) form a triangular lattice. Neutral sites 
with spin-| moments (Co 4+ ) form a Kagome lattice. Quan- 
tum fluctuations will reduce the amplitude of charge order 
substantially. Accompanying charge order, we expect some 
interesting spin liquid phase or complex short range spin or- 
der at low temperatures 

Before we consider influence of H2O we discuss some 
simple charge orders at x = j and | that are favored 
by electrostatics in the unhydrated Na x Co02 ■ We ig- 
nore the superexchange contribution, as J << Vi,V^. 
For x — j, the Co 3+ sites are arranged on a trian- 
gular lattice (figure 1) to minimize coulomb repulsions. 
Interestingly, the remaining sites carry spins and form 
spin-i Heisenberg antiferromagnet on a Kagome lattice. 
In our convention, the classical energy of this state is 
Ei (Kagome) = 0. In the real system, carrier dereal- 
ization will considerably reduce the amplitude of charge 
order. In this sense the charge ordered states shown in 
figure 1 and figure 2 are to be thought of as reference 
classical states. 

Another ground state comparable in energy is an 
anisotropic metal. It has ordered stripes - alternating 
insulating and 0.5 electron doped chains. The electro- 
static energy of this state per site is \{V\ + V2). How- 
ever, the carrier derealization in the 0.5 electron doped 
chains leads to a gain in kinetic energy which is easily 
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estimated when J is neglected in our t-J model. This 
case corresponds to a quarter filled infinite U Hubbard 
model, which can be converted into a half filled band of 
non-interacting spinless fermions giving us the dereal- 
ization energy = — |t|^cosfc = —2^-. Thus we get a 
total energy per site, Si (stripe) = j(Vi + V2) — 2^. 

Figure 2 shows the case of x = |. This classical ground 
state minimizes electrostatic repulsion and Co 3+ sites fill 
one of the three sublattices and the remaining hexagonal 
lattice is the neutral spin-| site. It is a hexagonal spin- 
5 quantum antifcrromagnet. The energy of this state per 
site is Ei = We also find striped states which arc 
local minima. 

x= 1/3 
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FIG. 2. 'Classical' charge order for x = |. Triangular lat- 
tice of localized charge — e carriers and a Hexagonal lattice of 
neutral spin-i moments. Quantum fluctuations will reduce 
the amplitude of charge order substantially. Accompanying 
charge order we expect short range AFM order at low tem- 
peratures 

So far we studied charge order in the C0O2 plane at 
commensurate fillings. As we move away from x = | 
and x = j, defects and discommensurations will be pro- 
duced or we may go to an incommensurate charge ordered 
structure. There may be one or more first order phase 
boundary between x = | and x = i. However, all these 
nice charge ordered phases will be challenged by a gener- 
ically disordered arrangement of Na ions in a triangular 
lattice and the consequent random potential seen by the 
mobile Co 3+ carriers as explained below . 

The energetically preferred sites of Na atoms [14] in 
Na x CoOi form a triangular lattice that have the same 
lattice parameter as the triangular Co layer. However, 
the Na and Co lattices are relatively shifted - if we 
project the allowed positions of the Na atoms of the 
nearest top and bottom layer onto the Co layer, these 
sites become the dual lattice (hexagonal lattice) of the 
Co triangular lattice. Because of this a sublattice order 
of the Na atoms does not couple to the charge density 
wave order-parameter of the C0O2 lattice and in principle 
allows a finite temperature charge order- disorder phase 
transition: a 3-state Potts (Z 3 symmetry) model transi- 
tion at x = ^ and a 4-state Potts model (Z 4 symmetry) 
transition at x = j. However, an inevitable disorder in 
Na sublattice leads to, based on an Imry-Ma type of ar- 
gument, a glassy order at low temperatures rather than 



a genuine charge order phase transition. Thus we ex- 
pect a phase diagram depicted in figure 3 for the range 
\<x<\. 

Let us move on to the hydrated case, N a x Co02-yH20 . 
As we mentioned earlier, the enhanced dielectric constant 
of the H2O layer will screen short range coulomb repul- 
sion and weaken and melt the high temperature charge 
ordering. It also screens and weakens the random Na po- 
tential seen by the carriers. For the appearance of low 
temperature superconductivity T c h need not be reduced 
to nearly zero value. A sufficiently weakened charge or- 
dered state may give up at low temperatures and super- 
conductivity may emerge. Figure 3 shows sketches the 
change of phase diagram as we go to the hydrated case. 
A strong resistance anomaly seen in a recent experiment 
[5] at T* « 50 K, may be the weakened charge order 
transition that we are discussing. 
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FIG. 3. Schematic phase diagram showing how hydration 
affects electronic phases: superconductivity is stabilized at 
low temperatures by a strong suppression of charge glass or- 
der. 



Let us discuss nature of H2O ordering hydrogen bond- 
ing in Na x Co02-yH20 in some detail. In what follows we 
use a recent suggestion of Cava et al. [15] that H2O may 
have a structure similar to the layers in hexagonal ice (lh 
ice). As mentioned earlier, in Na x Co02 , energetically 
favorable interlayer sites of Na form a triangular lattice. 
These sites are at the center of trigonal biprisms, capped 
by an oxygen atom at the top and one at the bottom. To 
understand H2O ordering, we consider Nai .Co02-^H 2 
as a reference compound. Fill one of of the three sublat- 
tices by Na atoms to minimize electrostatic energy. We 
are left with two empty sublattices that form a hexago- 
nal lattice. Two H2O molecules may be accommodated 
at the top and bottom of the capped trigonal biprism. 
By doing so we get two hexagonal lattices of H2O sand- 
wiching a triangular lattice of Na ions. Thus we have a 
triangular lattice filled by Na and H2O in the ratio 1 : 4. 
We can view the above hexagonal sheets as the sheets in 
hexagonal ice structure, as suggested by Cava and collab- 
orators [15] in their preliminary studies. We find it very 
interesting that the nearest neighbor H2O-H2O distance 
in the above geometry, w 2.81 Au, is nearly the same [16] 
as that in real hexagonal ice ~ 2.71. No wonder water 
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may freeze into ice in Na x C'o0 2 .yH 2 ! Having nearly 
the same H 2 -H 2 distance may also help ice sheets to 
have a good hydrogen bonding network like in hexagonal 
ice. 

Thus it is likely that H 2 molecules continue to have 
hydrogen bonding in spite of the Co0 2 and Na environ- 
ment. As hydrogen bond energy is substantial ~ 0.5 eV, 
water tends to have hydrogen bonding activity in ex- 
treme environments. Examples are biological systems, 
clathratc hydrates and water containing charged ions, 
where H 2 continue to maintain hydrogen bonding even 
though the local structure deviates considerably from the 
standard ice or water structure. Further the random 
Na environment in Na x Co0 2 .yH 2 may convert the 2D 
ice layer into a 2D amorphous ice layer with a good short 
range hexagonal order. 

Now we discuss how the dielectric property of the 
H 2 layer may control the low temperature electronic 
phases of the conducting Co0 2 layer. We are interested 
in finding how the electron-electron interactions at the 
charge order wave vector q = Q get screened by the 
interacting water dipoles. The relevant static dielectric 
constant is £h 2 o(Q)- ln an ice system like ours with a 
large disorder in dipole orientations, we expect very small 
variation of £h 2 o(<z) with q. Further random site disorder 
in the Na sublattice will produce Bjerrum defects; so we 
do not expect a 2D dipolar order-disorder transition, as 
in ideal 2D models of ice. 

Dielectric constant of ice has been studied extensively 
in the past and also recently [17,18,16]. A general expres- 
sion for the dielectric constant of an interacting dipolar 
system is [17,18]: 



+ 



47T 



Wk B T 



<(P-(P)) 2 



(3) 



Here ~ 1 — 2 is the high frequency dielectric constant 
of the dipole, in our case H 2 molecule. P is the total 
dipole moment of the system of volume V. And (...) 
denotes thermal average. Static dielectric constant of ice 
has not been measured at liquid He temperatures, as the 
dielectric relaxation becomes too slow even around liq- 
uid air temperatures. Fortunately, extensive numerical 
study of en 2 are available. For example, a recent calcu- 
lation [18] shows that for hexagonal ice, en 2 o ~ 220 at 
T = 50K. 

We use this 3D result to get an approximate estimate 
for our weakly coupled hexagonal ice layers as follows. 
We replace the volume V by « 6V to account for the 
c-axis expansion in Na x Co0 2 .yH 2 . Missing hydrogen 
bonds along the c-axis reduces the number of allowed pro- 
ton configurations leading to a reduction of ((P — (P)) 2 ) 
to w |((P — (P)) 2 ). This gives us en 2 o (hexagonal sheet) 
w 20 at T — 50K. Our system being strongly disordered, 
we do not expect en 2 o (hexagonal sheet) to change at 
lower temperatures. Thus the background dielectric con- 
stant of Na x Co0 2 .yH 2 is e = e a + e H2 o ~ 8 + 20. This 



reduces the mean short range repulsion by nearly a fac- 
tor of 3, making the T c h ~ IK, in equation (2). Once 
the long range charge order is disabled by a reduction of 
V, the simple t-J model and consequent low temperature 
superconducting phase is realized, albeit with a reduced 
T c in the range | < x < \. The sharp reduction in su- 
perconducting T c at the commensurate boundaries of the 
dome should arise from the strong short range order and 
lesser discommensurations and defects. The experimen- 
tally seen flat value of T c « 2K for x < \ and for x > |is 
likely to be an effect of phase separation. 
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FIG. 4. Schematic experiment to observe the 'Electrical 
Modulation of Superconducitivy'. The order voltage required 
is 500 V, when the film thickness is about a micron. 

Our proposal of a critical and catalytic role of 
H 2 layer suggests ways to access and control low tem- 
perature electronic phases of the Co0 2 layer. Based on 
this we suggest two effects: i) 'Electrical Modulation of 
Superconductivity'. Here we control (figure 4) super- 
conductivity and superconducting T c by modifying the 
screening property of H 2 layer by external electric fields 
- DC, AC or pulsed fields. This has interesting conse- 
quences of being able to locally erase superconductivity 
by STM tips, dynamically create Josephson networks or 
create 2D superconductivity of desired shapes through 
appropriate capacitor shapes etc. Since microwaves are 
absorved by hydrogen bonded networks we can pump 
microwaves at appropriate frequencies (hui < A sc , the 
superconducting gap) and dynamically polarize water 
dipoles and may influence its dielectric properties, and 
in turn control superconductivity. 

ii) 'Electroresistance Effect'. By modifying the am- 
plitude of charge glass order as well as T c h in the non- 
superconducting state, by influencing the H 2 layer by 
external DC or AC electric field, we can change p a b, the 
ab-plane resistivity. 

Below we estimate the electric field required to com- 
pletely suppress superconductivity. Having established 
a hydrogen bonded network it requires a finite energy 
to rotate a water molecule and orient its dipole moment 
along an external electric field. In infrared absorption 
and neutron scattering experiments [16] the absorption 
band corresponding to rotation of H 2 molecules is in 
the range 60 to 120 meV. Assuming a random orien- 
tation, the average energy required to reorient a water 
molecule is « 50 meV, i.e., a potential of 50 mV applied 
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over the length w 1 Au of the water dipole will orient the 
dipole moment of water along its field. If we have c-axis 
oriented Na x Co0 2 .yH 2 film of thickness 1 micron we 
need to apply a voltage w 500 Volts across the film, in 
order to orient the majority of dipoles. Strong polariza- 
tion of water dipoles reduces the dielectric constant, as is 
evident from equation (3). The resulting reduced screen- 
ing of carriers in the Co0 2 layer allows charge order to 
grow and superconductivity gets suppressed. 

A theoretical analysis, including some of the possible 
difficulties in observing the effects will be presented in a 
future publication. 

To get a clear understanding of this complex system, 
and to see if our proposal is correct more experiments are 
necessary: 

a) Charge Order: It will be interesting to perform 
NMR, NQR, STM, fiSR and other local probe measure- 
ments to look for charge order in the vicinity of the 
commensurate fillings x = 3 > g > 5 > § an d § and see how 
they differ between the two systems Na x Co0 2 .yH 2 and 
Na x Co0 2 ■ 

b) Spin order, singlets and gaps: Accompanying 
local charge order we expect a spin order at low tem- 
peratures (the scale of J is small, ~ 6 to 7 meV). In 
general the enhanced singlet stabilization by the superex- 
change process will introduce some kind of spin gap phe- 
nomenon. If the charge order at x = j leads to a Kagome 
lattice of spins it will be an interesting testing ground 
for some of the ideas of the spin liquid phase of spin- 
5 Kagome antiferromagnet, including possible novel ex- 
citations. 

c) Lower Doping: Experimentally, it has not been 
possible [3] to make Na x Co0 2 .yH 2 for x < \. It is 
likely [3] that c-axis ionic bonding is weakened, by de- 
creasing x and presence of water layer, making a 3D 
structure unstable. It will be important to synthesize, 
by non-equilibrium means, meta stable compounds for 
x < j to test the validity of RVB theory and also test 
our hypothesis of the role played by water. 

d) Replacing H 2 : It will be desirable to have a 
stable solid Na x Co0 2 -yX , where an intercalant 'X' not 
only increases the dielectric constant but also provides 
additional bonding between Co0 2 layers and make sta- 
ble compounds for x < \. 

e) Higher doping: According to reference 6, the 
dopant induced dynamics, within the t-J model will favor 
ferromagnetic correlations and a consequent p-wave su- 
perconductivity at higher dopings slightly above x = \. 
It will be interesting to look for this. 

f ) Inhomogeneous Superconductivity: If ice plays 
a central role, as suggested in this paper, H 2 density 
fluctuation in ice layer will directly influence supercon- 
ductivity in nearby Co0 2 layers resulting in a corre- 
sponding fluctuation in the superconducting order pa- 



rameter and possible well grown local charge ordered 
phase. 

g) Slow Relaxation: Since interacting water dipoles 
have very slow dielectric relaxation time scales [16], they 
may consequently affect superconductivity and impose 
some anomalous relaxation/aging effects. 

The present paper is phcnomcnological and qualita- 
tive in character. Any detailed quantitative calculations 
of T c and phase diagram for this complex system needs 
further experimental guidance. Issue of calculating lo- 
cal screening and dielectric constant in hydrogen bonded 
systems is known to have subtleties; added to this, we 
have conducting layers sandwiching water layers. We 
have made very crude estimates based on simple phys- 
ical arguments and very approximate considerations, as 
our primary aim is to focus and identify how water could 
play a crucial role in this complex system. 

We thank A.K. Mishra, V.N. Muthukumar, Debanand 
Sa, Manas Sardar and R. Shankar for discussion and 
Latha Malar Baskaran for reference [18]. 
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